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1  Introduction 


The  ability  of  high-resolution  direction  of  arrival  (1H)A)  estimators  such  as  Ml  SIC 
to  dist  inguish  between  closely  separated  sources  is  severely  degraded  bv  perl  ui  bat  ion 
errors  (the  term  pi  rl urbttl ion  error  refers  to  any  discrepancy  between  the  manifolds 
of  actual  and  theoretical  steering  vectors  for  the  array)  [li.  Some  possible  causes  of 
perturbation  error  are  element  misplacement .  uon-isot  ropicity  of  t  lie  elements  and/or 
the  environment,  inaccurate  gain  and  phase  calibrat  ion.  and  t  he  assumpt  ion  that  the 
sources  lie  in  the  far-field.  Various  methods  were  proposed  to  deal  with  tin*  problem 
of  estimating  DOA  under  perturbed  conditions.  Many  of  them  made  rather  restric¬ 
tive  assumptions:  [2]  required  at  least  two  calibrating  sources  with  known  DOAs. 
while  [3]-[  l]  required  t  he  sources  to  be  disjoint  in  eit  her  frequency  or  t  inir.  l  he  phase 
errors  were  assumed  to  be  caused  solely  by  element  mi-placement  in  [■">]-[(>].  while 
in  [T]-[9],  the*  phase  errors  were  assumed  to  be  independent  ol  angle  (and  therefore 
not  caused  by  element  misplacement).  1  his  report  describes  a  new  method  for  esti¬ 
mating  the  DOAs  of  noil-disjoint,  uncorrelated  narrowband  signals  which  makes  no 
assumptions  about  the  phase  errors  induced  by  the  perturbation  (like  other  methods 
which  include  phase  errors  in  tin  perturbation  model,  it  does  make  the  standard  as¬ 
sumption  that  the  gain  recurs  arc  angle-independent:  the  theoretical  gains,  however, 
may  be  non-isotropic  if  the  signals  are  far-field).  The  first  part  of  the  method  con¬ 
sists  of  an  iterative  algorithm  which  yields  estimates  of  tin'  signal  steeling  vectors. 
In  the  second  part  of  the  method,  these  steering  vectors  are  used  to  estimate  the 
DOAs.  This  two-part  structure  and  the  generality  of  the  perturbation  model  make 
the  new  method  similar  to  the  method  described  in  [It)]:  the  methods  have  not  yet 
been  compared. 

2  Estimating  the  Steering  Vectors 

Decomposing  the  signal  covariance  matrix 

The  algorithm  starts  by  obtaining  an  estimate  of  the  signal  covariance  matrix.  R.„. 
If  the  noise  covariance  matrix  is  assumed  to  be  a  multiple  of  the  identity  matrix,  an 
estimate  can  be  found  by  replacing  the  small  eigenvalues  id  t he  sample  data  covariance 
matrix  with  zeros  alter  subtracting  their  average  from  the  larger  eigenvalues 

Since  the  signals  are  uncorrelated.  the  total  output  signal  power  at  an  element 
ecpials  the  t  he  total  input  signal  power  (which  is  t  he  same  for  all  elements)  mult  iplied 
by  the  squared  gain  of  the  element.  Consequently,  the  gains  can  be  estimated  by  the 
square  root  of  the  diagonal  components  of  R„.  bet  ()t,  denote  the  estimated  gain  of 

the  A- 1 h  element,  for  k  =  1 . V.  Tin1  signal  covariance  matrix  is  gain-calibrated 

by  dividing  the  i,j th  component  of  the  matrix  by  ppp  for  all  i  and  j.  After  this 

For  ail  arbitrary  noise  covariance  matrix.  It,,  can  lie  estimate)!  a-  described  in  the  well-known 
paper  on  tin-  M  I  SI(  ’  a  Igor  it  lim  [  I  1] 


operation  has  been  performed,  all  elements  ean  be  assumed  to  have  unity  gain. 

Once  the  calibrated  Il„  has  been  obtained,  its  eigenvectoi  deeomposit  ion  is 
computed: 

R^aUSlU1.  (1) 

when? 

UfU  =  1.  (2) 

and  S  is  a  diagonal  matrix.  U  has  dimension  A  by  \l .  where  ,\l  is  the  number  of 
signals.  The  signal  covariance  matrix  has  another  decomposition. 

R„  =  vsv'.  (d> 

where  V  is  a  matrix  of  the  signal  steering  vectors  and  S  is  a  diagonal  matrix  ol  the 
signal  powers.  From  (1)  (d). 


where 


I  =  u  vs'7-'  s^-v'u. 


u  =  US-1/J. 


Since  UtVS1/2  is  a  square  matrix.  (  1)  implies  that 


I  =  s1/2vtu  utvs1/J. 


Since  S  is  diagonal. 


(UV,)  (U^Vj)  =  U  lot  1  <  i.j  <  .1/  and  /  ^  j.  (n) 

where  v,  denotes  the  / 1 h  column  of  V. 

Updating  the  Estimated  Steering  Vectors 

The  algorithm  attempts  to  find  a  set  of  vectors,  represented  by  the  columns  of  a 
matrix  V,  which  satisfy  (5).  V  can  be  initialized  by  setting  its  columns  ecptal  to  the 
theoretical  steering  rectors  that  correspond  to  peaks  in  Ml  SIC  (MUSIC  is  described 
in  more  del  ail  in  the  following  sect  ion ).  An  arbitrary  iterat  ion  ol  t  he  algorit  lint  begins 
by  re-estimating  the  first  steering  vector  while  the  other  estimated  steering  vectors 
are  kept  fixed.  1  he  first  transformed  steering  vector.  U  v(,  should  be  orthogonal  to 
all  the  other  transformed  steering  vectors.  I'lius.  it  should  lie  in  the  null  space  of 

t  rh  r  i  t  t  i  .1  .■  i  ii.1  .  i  i  .  i  i*  t  T  y  l  ■  l  ■ 


U  V],  wliere  V|  denotes  the  matrix  formed  by  deleting  the  lirst  column  of  V.  1  his 
null  space  has  a  dimension  of  !.  For  a  vector  z  in  the  null  space.  v(  should  satisfy 
the  under-determined  set  of  equations 

U'vi  =  oz  ((>) 


lor  some  arbitrary  scale  factor  n.  II  Vj  is  a  steering  vector,  it  should  lie  m  the  column 
space  of  U  (see  (1)  and  (:?)).  l  ints,  \q  should  also  satisfy 


v>  =  liq 


for  some  vector  cj.  Substituting  (7)  in  ((>)  and  making  use  ot  (2).  one  finds  that 


v,  =  oUS'/'z. 


{») 


After  V)  has  been  updated,  ot  her  steering  vectors  are  re-esl  miated  m  t  he  same  lashion; 
the  iteration  com  Inch's  with  the  update  ol  the  d/th  steering  vector  estimate.  It  is 
(dear  that  at  the  end  of  the  first  iteration  of  the  algorithm,  V  will  have  converged 
to  a  set  of  vectors  which  satisfy  (b)  and  lie-  in  the  column  space  of  U.  However,  as 
it  stands,  the  algorithm  is  incomplete.  Additional  structure  must  be  imposed  on  V. 
corresponding  to  the*  fact  that  tin*  absolute  value  of  each  component  in  V  is  assumed 
to  equal  1,  the  nominal  gain  of  flu*  calibrated  elements,  l  hus.  the  last  step  in  updating 
V)  consists  of  scaling  it;  its  phase  is  left  unchanged  but  its  magnitude  is  replaced  by  1. 
In  the  Z, 2-norm,  the  resulting  vector  is  the  one  with  the  correct  gains  that  is  nearest 
to  the  unsealed  v,.  After  the  scaling  operation,  vq  no  longer  satisfies  ((>)  or  (7). 
Consequently,  the  algorithm  must  be  applied  for  an  indefinite  number  of  iterations. 
Two  criteria  given  below  can  be  used  to  evaluate  the  convergence  of  the  algorithm; 
it  will  be  seen  that  the  algorithm  does  not  necessarily  converge  monotonieally  with 
respect  to  either  of  these  quantities.  The  average  orthogonality  of  the  transformed, 
scaled  vectors  is  defined  as 

Y.J2  (|/|(/(UtV,.  U*v/ ) 

'  T  J  —777 - 77  “  • 


where  the  angle  between  complex-valued  vectors  x  and  y  is  defined  as 


tui(j(x ,  y )  =  cos  1 


The  average  fractional  amount  that  the  scaled  vectors  lie  in  the  column  space  ol  U 
is  defined  as 

y-  |PU  V.l 

r  mi ' 

where  a  projection  matrix  Py  is  defined  by 

Px  =  X(XfXf ‘x1. 

The  algorithm  can  be  terminated  when  no  steering  vector  changes  significantly  from 
one  iteration  to  the  next,  or  after  a  fixed  number  of  iterations. 


3  Estimating  the  Directions  of  Arrival 

l'lie  theoretical  steering  vector  is  computed  under  the  assumption  of  no  gain  or 
phase  errors.  MUSIC  estimates  the  directions  ol  arrival  by  linelmg  the  values  ut 
0  which  yield  pe-aks  in  —10  h)g10(vj(l  -  P(  )v,;).  where  I  is  the  identity  matrix  of 
dimension  N  by  A.  The  corresponding  t  heore-t  ieal  steering  vectors  lie  tar  irom  the 
noise  subspace,  or  equivaie-ntly,  close  to  the  signal  subspace. 

Since  the  algorithm  described  in  the  first  section  provides  estimates  of  the  indi 
vicinal  signal  stee-ring  vectors,  one  can  independent  Iv  est  imate  l  lu  elirect  ion  of  arrival 
of  each  signal.  The  direct  itjn  ol  arrival  ol  the  j  t!i  signal  is  e-stimate-e!  by  line  ling  the 
single  value  of  0  which  maximizes  —10  lecglu(  %■*(  I  -  Pv  )v.;). 

4  Simulation  Results 

The  experiment  involved  a  circular  array  of  six  elements  with  A/2  spacing  between 
adjacent  elements.  The  noise  was  white  ami  its  variance  was  the  same  at  all  elements 
i.e.  the  noise  power  did  not  depend  on  th<-  element  gains.  The-  element  gains  wen' 
nominally  equal  to  unity  but  actually  varied  uniformly  from  0.5  to  1.5.  Both  MUSIC 
and  the  new  algorithm  used  the  gain-calibrated  R,s.  The  signal  din'd  ions  of  arrival 
were  0  =  — 10°,  0°,  and  10°;  each  signal  had  SNR  of  20  elB.  Twenty  trials  were  per¬ 
formed.  In  a  trial,  each  element  was  randomly  displaced:  tin*  x  and  y  displacements 
were  independent  and  uniform  over  ±.03A.  The  signal  steering  vectors  were  addition¬ 
ally  perturbed  by  angle-independent  phase  errors  of  up  to  .03A.  Once  tin*  perturbed 
steering  vectors  were  generated,  the  sample  covariance  matrix  was  randomly  formed 
from  500  snapshots.  The  theoretical  steering  vectors  were'  calculated  at  one-degree 
increments.  The  algorithm  was  allowed  to  run  for  200  iterations.  Only  once  in  the 
20  trials  did  MUSIC  have  three1  peaks  in  the  general  vicinity  of  the  true  direct ions- 
of-arrival.  The  mean  values  for  the  new  hum  hod's  direct  ion-of-arrival  estimate's  wore* 
—  1 1 .0°,  — 0.5°,  anel  10.5°.  The*  st  anelarel  eleviat  iems  wore'  1.8°.  2.1°.  a n el  1.5°. 

The  iollowing  figures  were'  ge'ne'iale’d  using  the*  eiata  bom  the'  twentieth  trial. 
Fig.  1  shows  that,  the  algorithm  improved  upon  the-  initial  MUSIC  estimates  of  the 
signal  steering  vectors;  each  estimate-el  ste-e-ring  ve-ctecr  was  eeanpaie-ei  wit  i i  the  cleisest 
actual  steering  vectevr  anel  t lie-  average  angular  se-paralion  was  determined  (the  gain 
errors  were  not  e-onsielereel  in  this  ealculat  ie>n ).  Fig.  2  shows  that  the-  t  rauslbrmeel 
estimated  stee-ring  vecte>rs  rapidly  be-remie'  e>rt  liogecnal  H>  eacii  eel  her  alter  a  few  itera¬ 
tions,  as  elesire-el.  This  is  one-  way  tec  e-stimate-  hear  we’ll  the-  algecrithm  is  performing. 
Tig.  3  reveals  a  se'conei  wav  to  estimate  the  algorithm  perloi  inane  e'  tlie  estimate-el 
steering  ve-eteers  should  lie  elose  tec  the  signal  subspaev.  Figs.  I  anel  5  illustrate*  the 
direction-ol-arrival  speet  ra  fe»r  MUSIC  and  the  new  met  heal,  in  t  he  e  ase  of  t  lie-  lntle»r, 
the  individual  curves  for  the-  three  signals  have-  beeni  superimpecse-d  ecu  the-  same'  pled  . 
Although  MUSIC  was  ue>t  capable  ecf  re-solving  the-  closely  space'el  signals,  tin-  three 
directiems-of-arrival  are-  apparent  in  lig.  5. 
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5  Summary  and  Observations 

This  report  lias  presented  a  two-part  met  hod  lot  est  uii.it  mg  t  lie  direet  ions  ol  arrival 
of  urn  orrelated  narrowband  sources  when  there  are  arbitrary  phase  errors  and  angle- 
independent  j^ain  errors.  1  he  signal  steering  vectors  are  estimated  m  the  first  part 
of  the  method:  in  the  second  part,  the  arrival  directions  are  estimated.  ]t  should  be 
noted  that  the  second  part  ol  the  method  can  be  tailored  to  incorporate  additional 
information  about  the  nature  ol  the  phase  errors,  lor  example,  il  the  phase  errors 
are  known  to  be  caused  solely  by  element  misplacement,  the  element  locations  can 
be  est  i  ma  t  ed  concurrent  I  v  with  the  1 )( ).  \  s  by  t  rying  t  o  match  tin-  I  heoret  nal  si  eeri  ng 
vectors  to  the  estimated  oiu's.  Smmiatioii  result-'  sudsiest  that,  for  genera!  pertur¬ 
bation.  the  method  can  resolve  closelv  spitted  sources  under  conditions  i  <  >  r  which  a 
standard  high-resolution  DOA  method  such  .is  Ml  S I ( '  fails. 
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